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Abstract 

Qh . We present a refined duality estimate for parabolic equations. This 

■^T estimate entails new results for systems of reaction-diffusion equations, in- 

cluding smoothness and exponential convergence towards equilibrium for 
equations with quadratic right-hand sides in two dimensions. For general 
C^ ' systems in any space dimension, we obtain smooth solutions of reaction- 

diffusion systems coming out of reversible chemistry under an assumption 
that the diffusion coefficients are sufficiently close one to another. 

> 1 Introduction 

O 



This paper presents a refined duality estimate for reaction-diffusion equations 



_i- ' arising in the context of reversible chemistry, of the form 



dtat - dt A x a t = (ft - a,) ( I Y[ a? - k J] af ] , i = l..n, (1) 



with the homogeneous Neumann boundary conditions 

V x a,i(t,x) ■ v(x) = for x G 50, t > (2) 

corresponding to the diffusion of n species Ai with concentration a, := a,i(t, x) > 
0, i = l..n at time t > and point x G 51 C K , each with its own diffusion 
coefficient dt > 0, and to the reversible chemical reaction 

ai^i + --- + a n A.^j8i^i + -"+j8nA. ) a„ft£N, (3) 
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where the above reaction is modelled according to the mass action law with 
the stoichiometric coefficients on, fti £ N and with the (constant) reaction rates 
I, k > 0. The mixture is assumed to be confined in a domain fl C M. N as implied 
by the homogeneous Neumann condition (2), where i/(x) denotes the outward 
normal vector to $7 at point x £ dfl. 

Moreover, wc systematically denote by £1t = [0, T] x ft (for any given T > 0) 
and by p' the Holder conjugate exponent of p, i.e. | + A = 1. 

The mathematical difficulties in proving existence, smoothness and large- 
time behaviour theories for systems like (1), (2) increase with the degree of the 
polynomials terms (appearing in the r.h.s. of (1)) as well as with the number n 
of equations or the space dimension N. 

The refined duality estimate that we shall derive in this paper depends yet 
on another parameter of (1), namely the maximal distance between the diffusion 
rates appearing in (1), that is 

S := sup {di} - inf {di}. (4) 

i=l..n t=l.. n 

It is easy to see that when 6 = 0, the system (1), (2) can be rewritten as the 
coupled system between n—1 heat equations (for the sums of two concentrations) 
and a single reaction-diffusion equation, which greatly simplifies the analysis 
compared to the general case when 8 > 0. In particular, the dynamics of system 
(1) for <5 = satisfies a maximum principle, which fails to be true for general 
diffusion systems with S > (except for special systems where the structure of 
the reaction terms enforces a maximum principle). 

Perturbation methods can sometimes be used in order to transfer at least 
partly the properties of the system with 5 = to the case when 6 > is small, 
see e.g. [HM]. 

Our new estimate is also particularly efficient in the case when S > is 
small, but it still gives results for all parameters S in some situations, and even 
when for example one of the diffusion rates is 0. Moreover, when the smallncss 
of d is required, it can be explicitly estimated. 

In order to obtain this estimate, we use a combination of ideas coming from 
maximal elliptic regularity, a Meyers-type estimate which provides an explicit 
perturbation argument, and duality methods in the line of e.g. [HMP, HM, 
PSch]. We end up with the following Proposition for parabolic equations with 
variable coefficients: 

Proposition 1.1. Let SI be a bounded domain ofM. N with smooth (e.g. C 2+a , 
a > 0) boundary dfl, T > 0, and p €]2, +oo[. We consider a coefficient function 
M := M(t,x) satisfying 

< a < M(t, x) <b< +oo for (t, x) e T , (5) 

for some < a < b < +oo, and an initial datum uq £ L P (D,). 



Then, any weak solution u of the parabolic system: 

dtU — A x (Mu) = on £It, 

< u(0,x) — uq(x) for iefl, (6) 

V x u-u(x) = on [O,T]x0fi, 

satisfies the estimate (where p' < 2 denotes the Holder conjugate exponent of p) 

h\\ L p(n T ) <(l + ^ a ,^)r 1/p ||«o||L P( o), (7) 

and where for any a, b > 0, q <E]1, 2[ 



Ca + b g 

D ^b,q '■= 7=, 2 ' b-a ' ( 8 ) 

'-,q — 



1 - c„+ b - 



provided that the following condition holds 



C^J-^Kl. (9) 

Here, the constant C m . q > is defined for m . > 0, q s]l, 2[ as f/ie fresi fiftoi is, 
smallest) constant in the parabolic regularity estimate 

W^xV\\li(CIt) < Cm.q II/IU^Ot): ( 10 ) 

where v : [0,T] x $7 — » R is £/ie solution of the backward heat equation with 
homogeneous Neumann boundary conditions: 

dfV + m A x v = f on CIt, 

v(T,x) = for xen, (11) 

V x u ■ v(x) = on [0,T]x<9n. 

We recall that one has C m _ q < oo for m > 0, q 6]1,2[ and in particular 
Cm. 2 < — (cf. Lemma 2.1 below). Note that the constant C m . q may depend 
(besides on m and q) also on the domain Q and the space dimension N , but 
does not depend on the time T. 

The consequences of Proposition 1.1 can be best understood in the case of 
the most standard reversible chemical reaction, that is when (3) writes 

A\ + A 3 ^ A2+A4 (12) 

and (1) becomes (after the rescaling of the nonessential constants k and I to 
unity) 

d t ai - di A x ai = (a 4 a 2 - 01 a 3 ), 

9 t a 2 - d 2 A x a 2 = -(a 4 a 2 - aia 3 ), 

d t a 3 - d 3 A x a 3 = (a 4 a 2 - 01 a 3 ), 

<9ta 4 — d 4 A ;c a 4 = — (04 a 2 — ai a 3 ). 



We recall that for this system (with the boundary conditions (2) and provided 
that di > for i = 1..4), existence of weak solutions in L 2 (logL) 2 was obtained 
in [DFPV] , together with the existence of strong (smooth) solutions when N = 
1. This result was improved by Th. Goudon and A. Vasseur in [GV] thanks to 
a careful use of De Giorgi's method, see e.g. [DeG]. They showed that strong 
solutions also exist when N = 2. We also refer to [CV], where smooth solutions 
were shown to exist in any dimension for systems with a nonlincarity of power 
law type which is strictly subquadratic, see also e.g. [Ama]. 

We also recall two results on exponential convergence to the equilibrium: 
First, exponential convergence in any H p norm in the one-dimensional case N = 
1 was obtained for the system (13) with boundary conditions (2) in [DF08]. This 
result is based on the use of the entropy/entropy dissipation method with slowly 
growing a priori L°°-bounds (cf. [TV], [DM]). It required "at most polynomially 
growing w.r.t. T" bounds for the quantities sup te r 0T i ||ai(t, Olli^fn)- A related, 
yet non-explicit approach to entropy methods for reaction-diffusion type systems 
can be found e.g. in [Gro, GGH]. 

In a later improvement [DFEqua], the authors showed exponential conver- 
gence to equilibrium in relative entropy avoiding any L°°-bounds on the solution. 
Thus, interpolating the weak global L 2 solutions constructed in [DFPV], one 
obtains exponential convergence towards equilibrium in any LP norm with p < 2 
for all space dimension N > 1. 

It is interesting to point out that for space dimensions N > 3 the existence 
of global classical solutions for general (even constant) diffusion coefficients and 
initial data is an open problem despite the fact that all L 2 solutions converge 
exponentially towards the constant equilibrium in LP with p < 2. Up to our 
knowledge, it is only known that if solutions to (13) with boundary conditions 
(2) would blow-up in the L°° norm, then such a concentration phenomenon 
would need to occur in at least two densities at at the same position xq 6 Q 
at the same time to > [HM]. Moreover, an upper bound on the Hausdorff 
dimension of singularities was given in [GV] . 

Thanks to Proposition 1.1, this paper is able to provide a direct proof of the 
result in [GV] (that is, without use of De Giorgi's method) when N — 2. More- 
over, we can obtain immediately the exponential convergence of the solution of 
(13), (2) towards equilibrium in L°° , which is a significant improvement on the 
above mentioned L 2_0 -convergence of [DFEqua]. It is remarkable that in this 
specific case, no smallness requirement for S appears in the assumptions. More 
precisely, we prove the 

Proposition 1.2. Let Vt be a bounded domain ofM 2 with smooth (e.g. C 2+a , 
a > 0) boundary dfl. For all [i = 1..4] assume positive diffusion coefficients 
di > and nonnegative initial data a^o £ L°°(Q). 

Then, there exists a weak nonnegative solution at £ L°°([0, +oo[xJ7) to the 
system (13) with homogenenous Neumann boundary conditions (2) subject to 
the initial data ato for all [i = 1..4]. 



Moreover, we denote for [i = 1..4] by aioo > the equilibrium values of the 
concentrations Oj." Thus, {<Xioo}i=i..4 is the unique vector of positive constants 
balancing the reaction rate 

Oloo a 3oo = 0-2oo Q4oo 

and satisfying the three (linear independent) mass- conservation laws 

aioo + a2oo = TqT / ( a io + a 2o) dx, 

aioo + aioo = TqT / ( a io + a 4o) dx, 

a2oo + a-ioo = foj" / ( a 20 + a 3o) dx. 
Then, there exist two constants K\,K2 > such that 

4 

V*>0, 5Z ll a *(*. O — «<oo||i~(n) <Kie" K2t - (14) 

»=i 

Moreover, the norm ||ffli||i,°°([o -foo[xfi) an d the constants Ki, K2 can be explic- 
itly bounded in terms of the domain fl, space dimension N, the norm ||fflio|[i°°(0) 
of the initial data and the diffusion coefficients di, [i = 1..4]. 



Proposition 1.1 also entails that a similar result to Proposition 1.2 holds in 
any space dimension provided that S > is small enough. More precisely, we 
prove the: 

Proposition 1.3. Let SI be a bounded domain ofM. N with smooth (e.g. C 2+a , 
a > 0) boundary dfl. For all [i = 1..4] assume positive diffusion coefficients 
di > such that < a = inf{di} l= i..4, < b — sup{di}j=i..4 and nonnegative 
initial data <Zjo G L°°(Q,) (i = \.A). 

Then, if 5 = b — a < 2 (C a+b 1+2 /n) > there exists a nonnegative weak 
solution a,i G L°°([0, +oo[xJ7) to the system (13) with homogeneous Neumann 
boundary conditions (2) subject to the initial data a^). 

Moreover (with the notation of the previous Proposition 1.2), there exist two 
constants K\, k 2 > such that 

4 

Vt>0, $]lk(*,0-aioo|U»(fi) <Kie~ K2t . (15) 

4=1 

Here, the norm \\ai\\L oo ([o,+oo[xn) an d the constants K\,K2 can be explicitly 
bounded in terms of the domain Q, space dimension N, the norm ||ajo||L°°(n) of 
the initial data, and the diffusion coefficients di, (i = 1..4). 



A further result based on Proposition 1.1 states an existence theorem for 
weak and (with more stringent assumption) bounded weak solutions of (1), (2) 
when the r.h.s. of (1) is not necessarily quadratic anymore. An assumption 
about the smallness of 6 > is still needed here. 

Proposition 1.4. Let fl be a bounded domain ofM. N with smooth (e.g. C 2+a , 
a > 0) boundary dCl. For all i = \..n assume positive diffusion coefficients 
di > such that < a = inf{dj}i=i.. Tl , < b = sup{di}»=i.. n and nonnegative 
initial data a^o S L°°(Q). Moreover, let k,l>0, cti,/3i € N be such that at least 
two coefficients fy — an are different from and have opposite signs. We define 
Q = sup{^™ =1 an, Yn=i @i} an d assume that Q > 3. 

Then, if S = b — a < 2 (Ca±b qi) (where Q' is the Holder conjugate of Q), 

there exists a nonnegative weak solution <Zj € L^(Q,x) for all T > to the 
system (1), (2) with the initial data a%Q. 

Moreover, if S = b — a < 2 (C a+b (o-i)(jv+2) )~ 1 , then the solution a, lies in 

2 > (Q-l) (JV+2J-2 

L°°{Q T ) for all T> 0. 

Finally, we show that there exist bounded weak solutions to the system (13) 
with homogeneous Neumann boundary conditions (2) in space dimension 2 even 
when one of the diffusion rates (say di w.l.o.g.) is equal to 0. We remark that 
existence of global weak solutions in the case of degenerate diffusion was also 
shown in [DFEqua]. 

Proposition 1.5. Let Q, be a bounded domain ofM 2 with smooth (e.g. C 2+a , 
a > 0) boundary dCl. For all [i = 1..4] assume nonnegative initial data e^o € 
L°°(Q) and nonnegative diffusion coefficients di > with e.g. di > 0, [i = 1..3] 
and e?4 = 0. 

Then, there exists a nonnegative weak solution ai lying in L°°(D,t) for all 
T > to the system (13), (2) (without the Neumann boundary condition on a^) 
subject to the initial data a^o- 



Remark 1.6. The assumption that the initial data lie in L°° enables to give a 
simple formulation of the Propositions above, but it is not optimal (if one is only 
interested in the bounds of the solutions after a given positive time tg > 0). For 
example, in the case of Proposition 1.2, it is easy to see that if the initial data 
lie in L P (Q) for some p > 2, then the conclusion remains true with the time 
interval [0, +oo[ changed into [to, +oo[ (for any io > 0) for the bounds. A more 
careful analysis (cf. Remark 2.3) shows that the assumption on the initial data 
can even be relaxed to L P (H,) for some p > 1. 

Remark 1.7. Classical bootstrap arguments also show that all the above weak 
solutions, once they belong to L°°(Qt), are in fact strong and smooth provided 
that the set O has a smooth enough boundary (and provided that the initial 



datum is also smooth enough, if one wishes to get smoothness even at point 
t = 0). Moreover in this case, those solutions are unique (in the set of smooth 
enough solutions) and (in the case of Propositions 1.2 and 1.3) converge towards 
equilibrium exponentially fast (with explicitly computable constants) in any H p 
norm (p G N), thanks to interpolation arguments similar to those exposed in 
[DF08]. 

Remark 1.8. The smoothness assumption C 2+a , a > on the boundary dfl is 
likely not optimal. One could conjecture that the above results hold true also 
for C 1+a boundaries or even Lipschitz boundaries. However, for more general 
boundaries, we lack, for instance, a reference which states explicitly the time- 
indcpendcncc of the constant C m . q in (10). 

The paper is organized as follows: In Section 2, we present the Proof of 
Proposition 1.1. Section 3 is devoted to the applications of Proposition 1.1 to 
the "four species" system (13), first in dimension 2 (Proof of Proposition 1.2) 
and then in any dimension (Proof of Proposition 1.3). Finally, we present in 
Section 4 the extensions to more general reaction-diffusion systems (Proof of 
Proposition 1.4) and to the case when one diffusion rate is (Proof of Proposi- 
tion 1.5). 

2 An estimate for singular parabolic problems 

We first recall a well-known result for the heat equation, which ensures that for 
m > 0, p G]l,2[ the constants C m , P stated in eq. (10) are well-defined, finite 
and time- independent. 

Lemma 2.1. Let il be a bounded domain of M. N with smooth (e.g. C 2+a , 
a > 0) boundary dfl, m > 0, and p e]l,2]. Then, there exists a constant 
C m ,p > depending on m, p, the domain Q, and the space dimension N , but not 
on T , such that the solution v : [0, T] x O — > R of the backward heat equation 
with homogeneous Neumann boundary conditions (11) satisfies 

||A x -y||j>(fi T ) < C„ hP \\f\\Lp(n T ), (16) 

where f € L p (VLt) is the r.h.s. of (11). Moreover, C TOi 2 < 1/m. 

Proof of Lemma 2.1. After introducing the time variable r = T — t £ [0, T] , the 
backward heat equation (11) with Neumann boundary conditions and zero end 
data transforms into the forward heat equation with zero initial data: 

d T v — m A x v = —f on Qt, 

v(0, x) = for x E fi, 

V x v ■ v(x) =0 on [0, T] x dil. 

Moreover, the semigroup of the forward heat equation with homogeneous Neu- 
mann boundary condition satisfies the contraction property ||e* x v(0, -)\\ p < 



||w(0,-)||p for all p £ [l,oo] and for all t > 0. Thus, the statement of the 
Lemma follows from [L], where it is explicitly stated that C m . p can be taken 
as time-independent. In particular, the Hilbert space case p = 2 allows explicit 
calculations by testing the above forward heat equation with —A x , which shows 
that C mt 2 < ^, see eq. (32) below. □ 

Wc now consider the corresponding problem with variable diffusion rate and 
obtain similar estimates: 

Lemma 2.2. Let ft be a bounded domain ofM. N with smooth (e.g. C 2+a , a > 0) 
boundary dfl. T > 0, p s]l,2], and M := M(t,x) be bounded above and below; 
i.-e. for some < a < b, 

< a < M(t, x) <b< +00 for (i, x) £ T - (17) 

We assume (using the notation of Lemma 2.1) that 

C^ ip t-1 < 1. (18) 

We consider f £ L p (Qt), and a solution v of the parabolic equation with 
variable diffusion rate given by M : 

dtv + M A x v = f on Q.t, 

< v(T, x)=0 for x £ O, (19) 

V x v-i/(x)=0 on [0,T]xdn. 



Then, 
and 



\\A x v\\ LP{nT) < D aAp \\f\\ LP{nT) , (20) 



IKO, Oilmen) < (l + ^-,6*)2 ,1 ^||/|| £ p ( n T ) ) (21) 

where D a ,6.p is given by (8), i.e. 



C< 



±b 



r\ ._ 2 -f 

^a,b,p •— b _ a , 

i-^4»,p— 

Proof of Lemma 2.2. In order to find estimates analogous to (16) for our vari- 
able coefficients parabolic equation, we take m := (a + 6)/2 and rewrite (19) as 
the perturbative problem 

d t v + m A x v = (m - M) A x v + /. (22) 

Then from (16), we get 

\\&xv\\ L p(n T ) <C m , P || (m- M)A x v + f\\ LP{nT) 

< C m , P ( h -^ \\^v\\ LP( n T ) + \\fhno T )) ■ ( 23 ) 



Provided that (18) holds, this directly implies (20). 
Using now estimate (20) in eq. (19), we get 

\\d t v\\ L p(n T ) < (l + &D a ,6j,)||/IUp(n T )- 
Taking into account that v(T, x) = for x G 0, 

rT 

d+v(t, ■) dt 

Lp{Q) 



(24) 



|v(0,-)IU-(n) 



< 



IMt, Oilmen) * < ||ft«(t, OIU^)^', (25) 



using Holder's inequality in the last step. Together with (24), this proves 
Lemma 2.2. □ 

From the previous Lemmas, we obtain by duality Proposition 1.1. 

Proof of Proposition 1.1. Take any / G LP (Qt)j and consider v the solution of 
the backward heat equation (19). Testing (19) with the solution u of (6), one 
easily checks that 



d 
dt 

which implies that 



u(t, x) v(t,x) dx 



u(t,x) f(t,x) dx, 



"f 



— / uq(x)v(0,x) dx 



< 



|wolU»>(n)IKMIIiy(n) 



< (I + bD aAp ,)T 1 /P \\u \\ LP{n) \\f\\ LP , {nT) , 



where we have used (21) for the last inequality, with p replaced by p' . As this 
holds for an arbitrary / G LP (CIt), we conclude that (7) holds, and Proposi- 
tion 1.1 is proven. □ 

Remark 2.3. In fact, one can observe that estimate (25) is not optimal. One 
can show using the properties of the heat equation that 

||«(0,-)IU'(O) < C T {\\d t v\\ LK n T) + \\A x v\\ LHUt )) , 

for any r < N f 9 _ 2 ; an d any T > [r can be taken as large as wanted if 
2p>JV + 2]. 

As a consequence, in the proof by duality of Proposition 1.1, the norm 
||uo||iP(fj) can be replaced by the weaker norm HwolU^n); f° r an y 9 > p/(l + 
2/N). This improvement allows one to consider more singular initial data in 
the reaction-diffusion problems studied in the sequel. 



3 The "four species" equation 

We now turn to the application of Proposition 1.1 to the "four species" sys- 
tem (13). 

3.1 A general a priori estimate 

We begin with the following a priori estimate for the "four species" equation, 
which is a direct consequence of Proposition 1.1: 

Lemma 3.1. Let VL be a bounded domain of R N with smooth (C 2+a , a > 
0) boundary dfl and T > 0. Consider a weak solution {fflt}i = i,...,4 to system 
(13) with homogeneous Neumann boundary conditions (2) on [0, T] and initial 
condition a^ G L p (il) (i = l.A) for some p > 2, and diffusion rates di > 
(i = l.A). We denote 



and assume that 



Then 



min {di}, b := max {di}, (26) 

i=l,...,4 i=l,...,4 



%v h -r < !■ ( 27 ) 



4 

IN| LP(nT) <(l + 6£> a , 6 . P 0r 1 /* 5>,o , (28) 

for i = l.A. 

Here, C TOjP and D a ^ iQ are the constants defined in Proposition 1.1. 

Proof of Lemma 3.1. We call u := Yli=i a i the total local mass of the system. 
Then u satisfies the forward heat equation (6) with 

Y" 1 - , di aAt, x) 
M(t,x):= ^ ;\ ■ (29) 

We observe that the bound (5) holds. Moreover, assumption (27) is identical to 
assumption (9). As a consequence, thanks to Proposition 1.1, we end up with 
estimate (28). □ 

Next we turn to a Lemma which is specially devised for the two-dimensional 
case. Note that it does not depend on the size of b — a: 



Lemma 3.2. Let Q be a bounded domain o/M 2 with smooth (C 2+a ) boundary 
d£l, T > and diffusion rates di > (i = l.A). We still use the notation (26). 



10 



Then, one can find two constants K\ > and p > 2 depending on a, b and 
Q, such that any weak solution (aj)j=i..4 of system (13) with the homogenenous 
Neumann boundary conditions (2) and initial conditions in L p (Jl) satisfies 



W^Wlv^^K^I + T) 1 ' 2 J2 a >o 



Lp(n) 



(30) 



Proof of Lemma 3.2. At first wc shall deduce that for any m > and for 3/2 < 
r < 2 the following csimtatc holds: 

C m ,r<m-$( r -i)(C m!3/2 )iV- r \ (31) 

Indeed, multiplying (11) by A x v and integrating on J7t, we easily obtain 

1 „ .„ 



||A a: w|| i 2 (OT) < 



|L 2 (n T ), 



(32) 



(33) 



so that C m> 2 <m l for m > 0. 

Then, an interpolation with (16) for p = | gives for all r <E [3/2, 2] 

\\A x v\\ L r {€lT) < m- e (C ra , 3/2 ) 1 - 9 ||/||L"(n T ), 

with the interpolation exponent e [0,1] satisfying 

9 1-6 1 4r - 6 

2 + 172" = r =* °-^-> 

which yields (31). 

Using now (31) for r = p' and m = ^rp, we see that eq. (27) is satisfied as 
soon as the following inequality is satisfied: 



-6(p') ln(m) + (1 - 6{p')) ln(C m , 3 / a ) + In 



b — a 



which yields with 1 - 0(p') = 4(2 -p') 

i(2-j/) (log^^+log^ 



<log 



<0, 



a + b 
b — a 



(34) 



Thus, the condition (27) is satisfied provided that ^^- Ca+b 3 , 2 > 1 as soon as 



3 logl^Cs*.^)' 



and, therefore, as soon as we choose ap'g [3/2, 2] satisfying 

1 l°gfe' 



2-f/< 



2 log(2±* C s 



L ,3/2 



Note that in the case s ^- Ca+b 3 , 2 < 1, the eq. (34) implies that condition (27) 
is always satisfied for all p' E [3/2, 2] since log(f^j) > 0. □ 



11 



3.2 Polynomial w.r.t. time bootstrap estimates 

We prove below a standard estimate for the heat equation (with Neumann 
boundary condition) which amounts to proving that the corresponding Green 
function has the same singularity as the Green function in the case of the whole 
a;-space R w 

We put the stress on the dependence of the constants w.r.t. the length T of 
the time interval (this dependence is only tracked with very great effort in the 
classical books like [LSU], cf. Remark 3.4 below). 

Lemma 3.3. Let f2 be a bounded domain ofM. N with smooth (C 2+a ) boundary 
dfl and T > 0. Let u be a solution of the forward heat equation with r.h.s. 
f £ L 9 (f^r) and homogeneous Neumann boundary conditions. 

dtu — d A x u = f on Qt, 

u(0,x) = uq(x) for x G ft, (35) 

Vu ■ v(i) = on [0,T]xdfl. 

If 1 < q < t 2 ' we consider s = < L +2 -2 ^ ^' an ^ assume further that the 
initial datum uq belong to L S (Q). 

Then, for any < e < 8 — 1, there is a constant Ct > depending on 
N, fl, e, d, q, |]uo||i s (o); ll/llL«(n T ) an d which has an at most polynomial 
dependence w.r.t. T , such that 

||tt|U—(n r ) < C T . (36) 

On the other hand, if q > i" , we assume that the initial datum uq belongs 
to L°°(Q). Then for any r G [1, +oo[, there exists a constant Ct depend- 
ing on N, Q, r, d, q, ||wolU°°(n); ||/||i«(fi T ) an d which has an at most 
polynomial dependence w.r.t. T, such that 



\Lr(n T ) 



< C T - (37) 



Remark 3.4. The statement of the above lemma is classical (and non even 
optimal) except that it crucially shows that the regularising effect of a parabolic 
equation involves constants which depend polynomially on the time interval 
[0,T] for all T > 0. In ID, this was already shown in [DF08] using a Fourier 
representation of the solution. For general domains however, the polynomial 
dependence of the constants seems to be nowhere in the literature. Moreover, 
tracking the constants, for instance, in the approach of [LSU] (where the Green 
function for the half-space problem are used locally along the sufficiently smooth 
boundary as transformed approximating problem) is much more difficult than 
the following proof. 

Proof of Lemma 3.3. We consider a solution of eq. (35). 

Step 1: Setting po = 9 an d testing (35) withpo |u| Po_1 sgn(u) (more precisely 
by testing with a smoothed version of the modulus \u\ and its derivative sgn(u) 
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and letting then the smoothing tend to zero) we obtain by integration-by-parts 
and Holders inequality with the constant Cq(po) '■= '■ 

!(ll<!°) +Co(po,d)Jjv x (u Po/2 )\ 2 dx<po\\f\\ L ?» (WuWH^ ■ (38) 



Then, the Gronwall lemma y <a(t) y ™ =► y{T) < y(0)»<> + ± J a(t)dt 
yields for all T > : 



i f T 



Ml 



NI%,CO< 



M| L ™ + 



) (s) ds 



< 2 P0_1 



hoir.c + n/ir.c^- 1 



< |_ll u o|| L ™ + ll/IL™ t fo 

'■= Ct,o- 



I'll 



We thus conclude that there is a constant Ct,o depending only onpo, d, \\uq\\ l p 
||/|| l po and polynomially on T such that 



sup \\u\\ p L ° P0 (t) < C T . . 
te[o,T] 



(39) 



Step 2: Gradient estimate and Sobolev embedding. The integration-in-timc 
of (38) and Holder's inequality show 



Cn 



v x K" /2 ; 



dx^WuoW^+poWfh.ojH^- 1 . 



The above estimate and Sobolev's embedding for H 1 with constant Cs yields for 
sq < oo for N = 2 and sq = jjh^, for iV > 2 (together with Young's inequality) 



T C 2 

\\u\\ p L \dt<^ 
Go 



pq- 1 



\\uo\\ p ; PO +Po\\f\\ L ?o(TC Tfi )— :=D Tfi , (40) 



where Dt,o is a constant depending only on p$, d, ||uo|Il po i II/IIl P0 an d P°ly~ 
nomially on T. 

In the Steps 3 and 4, we construct a sequence of exponents p n , s n and bounds 

(41) 

(42) 



sup 

te[o,T] 


\\u\\% r 


(*) < C T ,n, 


r T 








IMI& 


, dt < D T . n . 



In particular we set s n < oo if N = 2 and s„ = p w N _ 2 if JV > 3. 

Step 3: Iteration of (41): Similar to Step 1 we test (35) with p n +i u Pn+1 



-l . 



d 
~dt 



flN^Jii) +Cn+i / |v x K"+ 1 / 2 )| 2 dx=p„+i /" fu p ^- l dx, (43) 
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where C n +i(p n +i) := '■ p " +1 — — . In order to iterate the bound (41), we fix the 
exponent p n +i by introducing the n-independent exponent 

— -1 2 2 1 

r = — ^ : = 1 1 iff N > 3 which satisfy — < r < 1, (44) 

Sn - Pn+1 ^V Np p Q 

and any r satisfying — <r<lifiV = 2. Then, we estimate with p n +i — 1 = 
p n+ i(l — r) + s n (r — l/po) the above right-hand side of (43) by Holder's inequality 

f u ^(r-l/P0) uPn+ i(l-r) dx < iiju^ || u || S „(r-l/ P0 ) || u ||P 5 +^ 1 --) ) (45) 



and a Gronwall estimate for y < a(t)y r yields 



hW^tldT) < 



WuolCCXl+Pn+lT WfLroM^r^dt 



< 



|uo|| r P™+i +p„+ir||/|| L P C , / ||u 



,«»(r-l/P0)=^T 



dt 



l/r 



PI) - 1 



l/r 



Thus, by the definition of r we have s n (r — l/po) p< i 1 = s n -^ = p„, and we 
are able to use the bound (42) to obtain 



\«\\7&r(F)< 



ll"o||^ti+Pn+ir||/|| L ™^ 



1/7 



o 



T,n+1- 



(46) 



Step 4: Iteration of (42): Returning to (43) and (45), we collect 

2 



I 1( [IH| £ : + i]r ) + Cn+i ii tt iiS-^ +1 X h^) 



d.r 



<Pn+l ll/IL™ IMI L ^ 



»(r — - ) 



,p„+i(r-l) 

Holder's inequality as in (45) yield 



Since r < 1 we have ||tt|[ P i£+i 1; > (CT. n +i) r x and integration-in-time and 



{C T , n+ i) r - l rC n+l / V.K-+ 1 / 2 ) 



' dx<\\u4 r XX\+rpn + i\\f\\ L ™ D. p 



PQ- 1 
P 

T.i 



Finally, with s n +i = Pn+ijf^ ^ N > 2 and using Sobolev's embedding, 



Ikliptii * < 



C„+i r 



holl^+r^+ill/IL^i?^ 1 



D 



T.n+l- 



Step 5: Iteration in n. From the definition of r in (44) it follows that 
pn+i = s n (l — 7^-) + 7 and thus p n +i < 00 if N = 2 and for dimensions N >3, 
where s„ = Pnjjh} '■ 



N 



1 = Pr, 



J n+L — Fn 



M , 1 ^ (po - i) 



7V Po 



iV-2V W r "*po(JV-2)+2 p (VV-2) + 2 
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which has the fixed point Poa = jv+2 ^° 2 — and 

N + 2 N( Pa - 1) 

2 Po(jV - 2) + 2 

Thus, with pa ~ q and f <E L^ x we distinguish the cases 

J <7 < ^j^ where p„ — — 2h> p^ > q for <? > 1 , 
1 q > ^±2 where p„ "^°°> +oo. 



(47) 



In dimension N = 2 we can always choose p^ < +oo to be arbitrarily large. 
Note that for any n in the iteration, the constants Ct,u and Dt,u are polynomial 
with respect to T! 

Step 6: Interpolation of (41) and (42) in the cases poo < +oo (and thus 
N > 3). For any n we use Holder's inequality 

u p n u p n ir dxdt < [ |[u||£i„ ||u|| fp P „" (ft < C]! n D T , n . 
n T Jo 

In the limit n — > oo we find Poo^r = 7V+2-2« > 0- 

Thus, for all e > and in all dimensions N > 3 we obtain after finitely many 
iterations the following bound 

IMI P=0 «±2- e < Ct(IKIIl£=°,II/IIl? , q,d,c s ), 

t,x 

where Ct is a constant depending only on ||uo||£p°°, |[/||z,« , d, q, the Sobolcv 
constant Cs, and T. Moreover, Ct depends polynomially on T. □ 

Remark 3.5. We remark that 

1 1 N „, , 

^Eir-r;- 1+ —2 +0 ^ 

which corresponds to the regularity expected by convolution with the heat kernel 

N + 2 

being in L~ « ^, for all [i > 0. 

When applied to the quadratic "four species" eq. (13), the bootstrap above 
yields the following lemma: 

Lemma 3.6. Let £1 be a bounded smooth (C 2+a ) open subset ofM. N andT > 0. 
Consider then a weak solution {a.j},=i..4 to equation (13), (2) on [0,T] with 
initial condition {aio}i=i..4 € L°°(Q,t) and diffusion rates di > 0, [i = 1..4]. 

Assume that {aj}i=i..4 lie in L qo (p,T) for some <jo > (N + 2)/2, and that 
\\ a i\\ lio tn T ) grows at most polynomially w.r.t. T for i = 1..4. 

Then, for any r G [l,+oo[, we have {aj}j=i..4 G L t (Q,t) and ||ai||i7-(n T ) 
grows at most polynomially w.r.t. T. 
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Proof of Lemma 3.6. We use Lemma 3.3 repeatedly. In general, if a; G L q (Q,T) 
for some q > 2, then a.iOj G L q / 2 (Qt) {hi = 1, •••,4). Hence the right-hand 
side of eq. (13) is in L q / 2 , so from Lemma 3.3 we have 

f 1 9(.N+2) if1< </V + 2 

a t G L r - s (n T ) with r = 2 N + 2 ~i ^ 1< 9 < iV + 2, 

|oo ifg>A^ + 2, 

for any 6 > 0. If we define the sequence q n starting with the qo given in the 
Lemma, and satisfying 

9,1+1 = o m\ o _ — ' as long as 9n<N+ 2, (48) 

one can readily check that q n +i > <Zn is equivalent to q n > ^±^- and thus 
ln+l > gn+1 and wc obtain within finitely many iterations that 

qn+i q-n J J 

a, G L q "~ (Ot), with q„ > A^ + 2 for some n > 0, and any S > 0. 

Thus, applying once more Lemma 3.3, we end up with a,i G L t (VLt) for any 

r G [l,+oo[. □ 

In order to get an L°° estimate, we need one more computation: 

Lemma 3.7. Let 51 be a bounded domain ofM. N with smooth (C 2+a ) boundary 
d£l and T > 0. Consider a solution {aj}i=i..4 to the system (13), (2) on [0,T] 
with initial condition {aio}i=i..4 G L°°(Q,) and diffusion rates di > [i = 1..4]. 

Assume that {aj}j=i..4 is a weafc solution to the system (13), (2) on [0,T] 
satisfying a% G L 90 (J7t) /or z = 1..4 and some qo > (iV + 2)/2. ^4ko assume that 
|| a i||LTO(f2 T ) grows at most polynomially w.r.t. T. 

Then, 

IMU°°(fi T ) <Ct, i = 1..4, 

where Ct grows at most polynomially w.r.t. T . 

Proof of Lemma 3.1. Using Lemma 3.6, we know that dtai — di A x ai lies in 
L t (£It) for all r G [1, +oo[, with a norm in L t (D,t) which grows at most poly- 
nomially w.r.t. T. 

Using the result of [L], we obtain that the derivatives dtai and d XjXk ai lie 
in L t (Q,t) for all r G [l,+oo[, with a norm in L t (VLt) which grows at most 
polynomially w.r.t. T. 

A standard Sobolev inequality (in K. 7V+1 ) together with the use of an exten- 
sion/restriction operator implies that o^ lies in L°°(£It), with a norm in L°°(f2x) 
which grows at most polynomially w.r.t. T. □ 

3.3 Existence of bounded solutions and large time be- 
haviour for the "four species" model 

The results of the previous subsections enable us to prove Proposition 1.2 and 
Proposition 1.3. 
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Proof of Proposition 1.2. The existence of a weak solution to eq. (13), (2) [with 
given initial data in L°°(fl)] is already known (cf. [DFPV]). We observe that 
thanks to Lemma 3.2, we have a L p (Qt) estimate for a% (i = 1..4) for some 
p > 2. According to Lemma 3.7, the <Zj (i = 1..4) lie in fact in L°°(Qt)- Since T 
can be taken arbitrarily large, we end up with solutions defined on all R+ x it. 
Moreover, still according to Lemma 3.7, the L°°(£It) bounds of the a, (i = 
1..4) are at most polynomially increasing w.r.t. T. Using the cntropy/cntropy- 
dissipation estimate proved in [DFEqua] (or [DF08] , which used the assumption 
N = 1 only in order to show at most polynomially growing L°° bounds), we 
end up with the exponential decay towards equilibrium (14), and in particular 
we get a uniform- in-time bound for the a, in L°°(R + X it). □ 

Proof of Proposition 1.3. Once again, the existence of a weak solution to cq. 
(13), (2) [with given initial data in L°°(Cl)], is already known (cf. [DFPV]). 
Under the smallness assumption made on 5 = b — a, Lemma 3.1 implies an 
L p (Qt) estimate for en (i = 1..4) for some p > N/2 + 1. Then, Lemma 3.7 
ensures that the a% (i = 1..4) lie in fact in L°°(Ox)- The end of the proof (that 
is, the estimates about the convergence towards equilibrium) is exactly the same 
as in the Proof of Proposition 1.2. □ 

4 Extensions: General chemical kinetics and de- 
generate diffusion rates 

Proof of Proposition 1.4- We introduce the approximated system of equations 
constituted of the approximated equation (for any r G N*) 

(n n \ 
I FT (ap^ - k n (apft 
— ^i L, i = l..n, (49) 

i + H^-K) 2 ) 

where Q = sup{X]™ = i on, S™=i ft} ^ s defined as in the statement of the Propo- 
sition 1.4 and assumed to be supcrquadratic, i.e. we consider Q > 3. Here, 
cq. (49) is to be considered together with the homogeneous Neumann bound- 
ary condition (2) and a set of smooth approximated initial data {a[ }i = i : .. : „ 
(converging a.e. towards {dio}i=i,..,n as r — > oo and bounded in L°°(Q.)). 

The existence of a smooth (strong) solution to this approximated system 
follows from standard existence results of systems of reaction-diffusion equations 
with bounded and Lipschitz-continuous r.h.s. (cf. [D, QS], for example). 

With the assumption that at least two coefficients on — fy are different from 
zero and have opposite signs, one can find coefficients 7^ > such that 



£>(a*-ft) = 0. 



(=1 
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At first, we observe then that 

n 

9 t (£>al)-A x (M r <) = 0, 

where 

E i= i 7* <(*,#) 
is bounded a < M r (t,x) < b a.c. in £It- 

We now assume that 

b-a<2(C^ Qt )-\ ^ + ^7=1- (50) 

Using Proposition 1.1, we see that for any i — 1, ..,n, the sequence (a£) r gN i s 
bounded in L®(Qt) (f° r au T 1 > 0). Moreover, since (50) is a strict inequal- 
ity, an interpolation argument similar to the one used in Lemma 3.2 implies 
that (a[) rS N is bounded in L® +6 (Qt) for some sufficiently small e > 0. As a 
consequence, the quantities <9 t a[ — di /S. x a\ are bounded in L 1+£ 'Q{Q,t). 

The sequence (a[) r6 N converges therefore (up to extraction of a subsequence) 
a.e. as well as strongly in L®(£It) towards a limit ai E L^ +e {Qt) ■ Finally, one 
can pass to the limit r — > oo without difficulties, which ensures that the limiting 
concentrations a, satisfy the original system in the weak sense. 

Secondly, if we suppose 

— CI < 2 (C a + b (Q-l) (JV + 2) ) , 

2 I (Q-l) (JV + 2)-2 

Proposition 1.1 (and the same interpolation argument as previously) ensures 
that the weak solution defined above satisfies the extra estimate a% G L z °(Qt), 
for some z > (l+N/2) (Q-l). Then, d t ai-diA x ai E L z ^ Q (fl T ), and thanks 
to the properties of the heat kernel summarised in Lemma 3.3, a, G L p {VLt) for 
all p < zi, with z\ = „ z % (or all p E [1, +oo[ if z Q > Q (1 + N/2)). 

^ N + 2 Z(i 

We see that a simple bootstrap ensures that a% E L p (£It) for all p < Zk, with 

z k = n Z %~ 1 ( or an P e [Ij +°°[ if z k~i > Q (1 + N/2)). The sequence 2^ 

is increasing up to the point when Zk > Q (1 + N/2), therefore we obtain the 
estimate a* G L p (£It) for all p G [1, +oo[. We proceed as in Lemma 3.7 to get 
the final estimate a* G L°°(£It)- d 

Proof of Proposition 1.5. Existence of weak solutions (in L 2 (Q,t)) for the set of 
equations considered in this Proposition is shown in [DF07]. 

By adding the equations satisfied by a\ and a 2 , we see that 

d t (a 1 +a 2 )-A x (M(a 1 +a 2 )) = 0, (51) 

where M(t, x) E [inf{di,d2},sup{<ii, d 2 }] almost everywhere. 
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Adding the equations satisfied by a 2 and a 3 , we also see that 

d t (a 2 + as) - A X (M (a 2 + 03)) - 0, (52) 

where M(t, x) G [inf{<i2,<i3},sup{<i2, ^3}] almost everywhere. 

As a consequence, we see that thanks to Proposition 1.1 and an interpolation 
argument similar to the one used in Lemma 3.2, for some S G]0, 2[ (and any 
T > 0), a z G L 2+s (n T ), when i = 1,2, 3. Then, 01 a 3 G L 1+s / 2 (n T ). Since 

d t a 2 - d 2 A x a 2 < a±a 3 , 

we see using the properties of the heat kernel in 2D (cf. Lemma 3.3) that 

2 + S q 

a 2 G L 1 - 6 ' 2 (fir) (here and in the sequel, the notation L p ~ v means n 9e [i jP [L 9 ). 

Using a duality estimate (see e.g. [PiclO, Lemma 3.4], [QS, Lemma 33.3]) 
for solutions of (51), it follows that if a 2 G L 9 (£7t) for any 1 < q < 00 then 
also a\ lies in L q (JlT)- Thus, we deduce from the estimate on a 2 that 01, 

2 + S q 

a 3 G L 1 - 6 / 2 (CIt)- Wc build the (finite) increasing sequence p n G]2,4[ such 

that pn = 2 + S, and ^— = — I. We denote by Nn the last index such 

that pn < 4. The properties of the heat kernel in 2D (once again and in all 
the sequel, cf. Lemma 3.3 for a precise exposition of those properties) and the 
duality estimate of [QS] implies that <n 6 L PN o+ 1 (Q T ) ] when i = 1,2,3. A last 
application of the properties of the heat kernel in 2D and the duality estimate 
of [PiclO, QS] shows that a, G L°°~°(£It), when i = 1,2, 3. Finally, thanks to a 
computation similar to the one in Lemma 3.7, a 2 G L°°(51t). 

Observing that 

d t a,4 < aia$, (53) 

we also see (performing the integration in time) that 04 G L°°~ (f2T)- Then, 
for i = 1,3, 

d t a.i ~ diA x a,i < a 2 a 4 , 

so that thanks to the same computation as above (similar to the one in Lemma 3.7), 

ai,a 3 Gi°°(r2 T ). 
Using once again eq. (53), we obtain that 04 G L°°(VLt)- □ 
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